Mon. Not. R. Astron. Soc. 000, [TV?? (2011) Printed 9 August 2012 (MN WI&L style file v2.2) 



Evolution of the Population of Very Strong Mg II Absorbers 

. - . Paola Rodriguez Hidalgo 1 ' 2 * Kaylan Wessels 1 ' 3 , Jane C. Charlton 1 , 
y—i ■ Anand Narayanan 1 ' 4 , Andrew Mshar 1 , Antonino Cucchiara 1 ' 5 , and Therese Jones 1 ' 6 

■ ^ . 1 Department of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802 

V N ■ 2 Department of Physics & Astronomy, York University, Toronto, ON, Canada M3J IPS 

{^JQj 3 Materials Department, University of California, Santa Barbara, CA 93106 

r—{ _ 4 Department of Earth & Space Sciences, Indian Institute of Space Science & Technology, Thiruvananthapuram, India 695547 

, ^Department of Astronomy and Astrophysics & UCO/Lick Observatory, University of California, 1156 High Street, Santa Cruz, CA 

■ 95064 

QQ 6 Department of Astronomy, University of California, B-20 Hearst Field Annex, Berkeley, CA 94720 



Q • 9 August 2012 

u; 

43 ABSTRACT 
Oh 

q . We present a study of the evolution of several classes of Mg n absorbers, and 

$h 1 their corresponding Fe n absorption, over a large fraction of cosmic history: 2.3 to 

! 8.7 Gyrs from the Big Bang. Our sample consists of 87 strong (W r (Mg n) > 0.3 

■ A) Mg II absorbers, with redshifts 0.2 < z < 2.5, measured in 81 quasar spectra 

obtained from the Very Large Telescope (VLT) / Ultraviolet and Visual Echelle Spec- 
trograph (UVES) archives of high-resolution spectra (R ~ 45,000). No evolution- 
ary trend in W r (Fe n)/W r (Mg n) is found for moderately strong Mg II absorbers 

Q\ . (0.3 < W r (Mg n) < 1.0 A). However, at lower redshifts we find an absence of very 

■ strong Mg n absorbers (those with W r (Mg n) > 1 A) with small ratios of equivalent 

widths of Fe II to Mg II. At high redshifts, very strong Mg II absorbers with both 
small and large W r (Fe n)/W r (Mg n) values are present. We compare our findings to 
a sample of 100 weak Mg n absorbers (W r (Mg n) < 0.3 A) found in the same quasar 
spectra by Narayanan et al. (2007). 

The main effect driving the evolution of very strong Mg n systems is the difference 
between the kinematic profiles at low and high redshift. At high redshift, we observe 
that, among the very strong Mg II absorbers, all of the systems with small ratios of 
W r (Fe n)/W r (Mg n) have relatively large velocity spreads, resulting in less saturated 
profiles. At low redshift, such kinematically spread systems are absent, and both Fe II 
and Mg n are saturated, leading to W r (Fe ll)/W r (Mg n) values that are all close to 
1. The high redshift, small W r (Fe n)/W r (Mg n) systems could correspond to sub- 
DLA systems, many of which have large velocity spreads and are possibly linked to 
supcrwinds in star forming galaxies. In addition to the change in saturation due to 
kinematic evolution, the smaller W r (Fe Il)/W r (Mg n) values could be due to a lower 
abundance of Fe at high redshifts, which would indicate relatively early stages of star 
formation in those environments. 

Key words: galaxies: evolution — halo — intergalactic medium — quasars: absorp- 
tion lines. 



1 INTRODUCTION 

Quasar absorption lines allow the study of galaxies and their 
E-mail: prh@yorku.ca halos with no bias toward specific environments because of 



(N 



X 



© 2011 RAS 



2 P. Rodriguez Hidalgo et al. 



the random distribution of lines of sight in the sky. In partic- 
ular, high resolution studies of quasar absorption lines help 
us better characterize the kinematic properties and chem- 
ical content of the absorbing gas. The kinematic composi- 
tion of individual galaxie s can be studied statisticall y with 
quasar absorption lines (Ch arlton fe Churchilll [l998l Pro- 
file shapes of the absorption lines show clouds of gas at dif- 
ferent velocities relative to galaxy centroids, establishing, 
for example, whether the absorption is consis tent with ab- 
sorption in the disk and /or the extended halo l|Steidel et al.l 
l2002l ; lKacprzak et alll20l"lh . Also, understanding the nature 
and evolution of galaxies requires assessing their star for- 
mation h istory, which ca n be traced by their metal con- 
tent fe.g.. lMatteucci|[2008l and references therein). Although 
the number of massive stars is smaller than the number 
of less massive stars, the former have a larger impact on 
the galaxy's chemical enrichment through supernovae ex- 
plosions. Different types of supernovae (SNe) contribute dif- 
ferently to this enrichment: core-collapse SNe (Type II) pro- 
duce a larger contribution of a-elements (such as Mg) rela- 
tive to Fe within the first billion years of the formation of 
a stellar population. Fe is mostly generated in a later phase 
by Type la SNe. Thus, the ratio of Fe relative to Mg may 
be used as a measure of the stage in the star formation his- 
tory of a galaxy. Quasar absorption line studies are a probe 
of this kind of information. 

Among the absorption features that can be used in 
extragalactic studies, the resonance doublet transitions of 
one a-element, Mg n AA2796,2803, are well suited due to 
their strength and their accessibility at optical wavelengths 
from redshifts of 0.2 up to 2.7, covering a large range in 
cosmic ages (from approx. 2.5 to 11.3 Gyrs). These inter- 
vening Mg II absorption lines have been historically (and 
originally arbitrarily) classified based on the strength of 
their rest-frame equivalent width in the 2796 A transition 
(Wr(Mg Iij3). Posteriorly, physical properties and origins 
for these systems have defined statistically distinct popu- 
lations, although the equivalent width boundaries between 
these different populations are not sharp. Strong Mg II ab- 
sorption systems are those with W r (Mg n) > 0.3 A, while 
weak Mg II absorption is defined to have W r {Mg n) < 0.3 
A. Several sub-classifications have been proposed among 
the strong absorption systems: very strong absorbers are 
those with W / r (Mg n) > 1 A, and the term "ultra- 
strong" is reserved f or those Mg H absorption systems wit h 
W r (Mg n) > 3 A (|Nestor et al.ll2005l : INestor et al]|2007f ). 
Strong Mg II absorption has generally been found to be 
connected to galaxies. Luminous galaxies (~ 0.1 — 5L* 
galaxy) near quasar lines of sight show st rong Mg H ab- 
et al.l 
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halos are patchy, with ^50% covering fraction for Mg H 
llTripp fe Bowenll2005l : IChurchill et~aill2007l ; iKacprzak et all 
2008), and it seems this covering fraction might be even 
smaller at smaller re dshifts (less than 40% f or strong Mg II 
absorbers atz~ Q.l; |Barton fe Cookelfeoogl '). 

Strong Mg II absorbers have been found t o be good 
tracer s of cold/warm and low-ionization gas. iRao et all 
(2006) suggests that » 80% of very strong Mg II absorption 
traces damped Lyman-a (DLA) galaxies (JV(H i) > 10 20 ' 3 ), 
but the relation between A r (H i) and iV(Mg n) doe s not 
prevail for the largest W r (Mg ill INestor et all (|2007T l indi- 
cates that "while it is likely that a large fraction of ultra- 
strong Mgll absorption (perhaps > 60%) have iV(H i) > 
10 20 - 3 atoms cm , most (perhaps ~ 90%) DLAs have 
W r (Mg n) < 3 A". Indeed, W r (Mg n) and N(¥L i) do not 
correlate as much in this regime and ultrastrong Mg II sys- 
tems can be found in sub-DL A galaxies (10 19 < A r (H i) < 
10 20 3 atoms cm" 2 ) as well (|Rao et all l201lt ). Moreover 



1 Throughout the paper we use W r (Mg II ) to describe the rest- 
frame equivalent width of the Mg II A2796 transition 



iKulkarni et"afl (|2010h has found that the Mg II associated 
with sub-DLAs tends to show larger velocity spreads on 
average than that associated with DLAs. Together with 
the finding tha t sub-D LAs are more metal-rich than DLAs, 
IKulkarni et"aTl (|2010h suggests that the difference is due to 
their star formation histories, where galaxies associated with 
sub-DLAs would tend to be more massive, and suffer faster 
star formation and gas consumption, which would result in 
less N(H i) and more metal production. If the star forma- 
tion has occurred recently, we would expect to see kinematic 
signatures of it in the Mg n absorption profiles, due to out- 
flows, also called superbubbles or superwinds. Moreover, we 
would expect evolution of very strong Mg II absorption pro- 
files as the star formation rate in galaxies shows a decline at 
redshifts lower than z ~ 1. 

Indeed, Mg II absorbers have been observed to show 
evolution with redshift. In a survey of ~ 3700 quasars from 
an early dat a relea se of Sloan Digital Sky Survey (SDSS), 
INestor et all l|2005h found that the number of very strong 
Mg II absorbers (W r > 1 A) was not consistent with the 
expectations for cosmological non-evolution, showing larger 
numbers at high redshift (z > 1 .2). This trend i s mos t sig- 
nificant as W r (Mg n) increases. iProchter et all (|2006t ) con- 
firmed these results for the incidence of Mg II in a study of 
45,023 quasars from the Data Release 3 (DR3) of SDSS. The 
decline in number of very strong Mg II absorbe rs with time is 
consi stent with the decline of cosmic SFR fe.g. jMadau et all 
1998), which is found to have decrea sed by about one or- 
der o f magnitude since its peak (e.g. lRosa-Gonzalez et all 
|2002| and references therein). Since very strong absorbers 
are evolving away at redshifts coincident with the decrease in 
SFR, they could be tracing th e star formation sites of galax- 
ies, as previously suggested bv lGuillemin fe Bergeron! (Il997l ) 
and IChurchill et all (|l999h . Low -m ass galaxies present th is 
peak of SFR at even lower redshifts (|Kauffmann et al]|2004f ). 
In a study of z < 0.8 Mg 11 absorption of SDSS DR3, 
iBonche et all f2006) found an anti-correlation between halo 
mass of galaxies and the equivalent width of Mg 11 absorp- 
tion present. Since the very strong Mg 11 systems may show 
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saturation in most of their profiles, the equivalent width is 
determined primarily by the velocity dispersion of the ab- 
sorbing clouds. The required velocity spreads are consistent 
with a starburst picture for the strongest Mg II systems, 
but not with s tructu re within individual virialized halos. 
IProchter et~ai1 l|2006l ) supported the claim that very strong 
Mg II absorbing structures are related to superwinds rather 
than accreting gas in galaxy halos, in a study over a larger 
redshift range (0.35 < z < 2.7) and based on the kinemat- 
ics of W r (Mg n) > 1 A absorbers. High resolution studies 
of individual absorbers confirm that they show particular 
signatures in their profiles that might be consist ent with su- 
perwinds (|Bond et al.ll200ll ; lEuison et al.ll2003r i. Moreover, 
field imaging of a subset of the strongest Mg II absorbers 
(W r (Mg n) > 2.7 A) at low redshift (0.42 < z < 0.84) 
indicates that interactions, pairs, and starburst rela ted phe- 
nomena are likely to be present (jNestor et al.ll2007h . 

Weak M g n absorption systems hav e been found to 
evolve as well INaravanan et al.l ()2005l ) and lNaravanan et al.l 
found a peak at z = f.2 in the number per unit 
redshift, dN/ dz, of Mg u absorbers w ith W r (Mg n) < 0.3 
A. Moreover, INaravanan et akl l|2008l ) find a trend in the 
ratio of Fe II to Mg II: at higher redshift (z > 1.2), 
weak Mg II absorption systems with large values of 
W r (Fe llA2383)/VK r (Mg n) are rare. They suggested that 
this trend could either be caused by an absence of high den- 
sity, low ionization gas at high-z or the absence of enrich- 
ment by Type la supernovae in weak Mg II clouds at high-z. 
They suggest that the relatively few weak Mg II absorbers 
that are observed at high z are from young stellar popula- 
tions and thus a-enhanced. 

This paper investigates the evolution of very strong 
Mg II absorption systems with redshift, particularly their 
Fe II to Mg II ratios. We enlarge the s a mple of strong Mg II 
absorption presented by iMshar et alj (|2007 ), and compare 
it to t he analysis previously carried out in INaravanan et al.l 
l|2007t ) for weak Mg II absorbers. The data and methods, as 
well as the systems found, are presented in section ij2] Re- 
sults involving evolution of the profiles and the Fe n/Mg II 
ratio are shown in section §3] We summarize the results and 
discuss their implications in section ^4] 



2 DATA AND SURVEY METHODS 

We analyzed 81 high resolution R ~ 45,000 (= 6.7 km s -1 ) 
VLT/UVES quasar spectra retrieved from the ESO archive. 
This particular set of spectra was first compiled for the 
INaravanan et al.l <|2007t) search for weak Mg n AA2796, 2803 
absorbers, although the spectra were originally obtained to 
facilitate a heterogenous range of studies of the Lya forest 
and of strong metal-line absorbers. Because some lines of 
sight were observed because of previously known strong ab- 
sorbers, we expect a bias toward strong Mg II absorption. 
Since our focus is on the evolution of absorber properties, 
and not on tabulating a statistical distribution function of 
equivalent widths, this bias does not affect our study. How- 



ever, we remain alert to the fact that a particular kind of 
strong absorber could be favored in this study if commonly 
observed by the original studies. The list of quasars used, 
with inform ation on each qua s ar obs ervation, can be found 
in Table 1 of Narayanan et al. ll2007h . The data reduction is 
described in INaravanan et al 120071 . section 2.1. The average 



signal-to-noise ratio of the spectra is large (> 30 pixel -1 over 
the full wavelength coverage), which facilitates detection of 
weak absorption components. 

Although UVES spectra have a large wavelength cov- 
erage (3000 A to 1 /im), the choice of cross-disperser set- 
tings may cause breaks in this coverage. We searched the 
same r edshift path lengt h avail able in each quasar spectrum 
as did INaravanan et all (|2007T ) (illustrated in their Figure 
1), which excludes the Lya forest and wavelength regions 
contaminated by telluric features. We also exclude systems 
within 3,000 km s -1 of the quasar redshift in order to avoid 
contamination with intrinsic associated absorption lines. Al- 
though strong systems are easier to spot than weak systems 
in contaminated regions, we avoided such regions because 
we could easily miss weaker subsystems associated with such 
systems, which would bias our resultfl In fact, strong ab- 
sorption systems are almost always composed of more than 
one subsystem. 

We searched for strong Mg 11 (W r > 0.3 A) systems 
along the 81 quasar lines of sight and found a total of 87 
systems along 47 of the sight lines. We used a 5a search 
criterion for Mg 11 A2796 and at least 2.5a for Mg 11 A2803, 
which was sufficient for strong Mg 11 absorbers, especially 
due to the high S/N of the spectra. After confirming the 
presence of Mg 11 AA2796, 2803 by visual inspection of the 
profile shapes, we looked for other ions that could be present 
at the same redshift: Fe 11 A2374, A2383, A2587 and A2600, 
and Mg I A2853. The Fe 11 and Mg I transitions can be 
used for a better understanding of the kinematics of sys- 
tems that appeared saturated in Mg II. Since these ions are 
weaker than Mg 11 we searched for them using a 3a detec- 
tion limit. We also compared the profile shapes of the Fe II 
transitions and of Mg I to those of Mg 11 AA2796, 2803, in 
order to identify possible blends contaminating the profiles. 
We placed upper limits on equivalent widths whenever such 
blends occurred. Equivalent widths of all detected transi- 
tions were computed by a pixel by pixel integration of the 
profiles, including all subsystems, as in CV01. The absorp- 
tion redshift of each system was formally defined by the 
optical-depth-weighted median of the Mg 11 A2796 profile, 
following Appendix A in CV01. 

The 87 strong Mg 11 systems are shown in Figure Set 
[1] and listed in Table [1] Figure Set [1] includes profiles of 
Mg 11 AA2796, 2803, and of Mg 1 A2853, Fe 11 A2374, Fe 11 



2 We define 'subsystem', following the CV01 definition, as 'ab- 
sorption features that are separated by more than 3 pixels (i.e., 
a resolution element) of continuum flux'. Mg II doublets found 
within 1,000 km s — 1 of each other were classified as subsystems 
that are part of the same absorbing system. 
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A2383, Fe n A2587 and Fe n A2600, whenever the spectral 
region was covered. In these figures, the spectra are nor- 
malized and the different transitions are aligned in velocity 
space, centered at the median of the apparent optical depth 
of the Mg n A2796 profile. We included vertical dashed lines 
to delineate the separate subsystems, although our equiva- 
lent width measurements are only tabulated for the systems 
as a whole. 

Table [T] lists the rest-frame equivalent widths or equiv- 
alent width limits for the Mg II, Fe II, and Mg I transitions 
of the 87 strong systems. Among them, 58 systems are mod- 
erately strong (0.3 A < W r (Mg n) < f A) and 29 systems 
are very strong absorbers (W r (Mg n) > I A). Blanks in Ta- 
ble Q] represent cases where the wavelength of the relevant 
transition was not covered in the VLT/UVES spectrum. In 
cases where the Fe II transition or Mg I was covered, but not 
detected, 3<r upper limits are given. We visually inspected 
every Fe II absorption feature, and carefully compared their 
profiles to each other and to the Mg II profiles, in order to 
confirm that they were consistent. We rejected deviant parts 
of these components, mostly caused by blends with different 
transitions from systems at other redshifts. In cases where 
blends prevented the measurement of an accurate equiva- 
lent width, we place lower limits by avoiding the blended 
region and upper limits by including it. These results ap- 
pear as a range of values in Table [T] (e.g., 0.354-0.426 A 
for the Wi-(2587) measurement for the absorption system at 
z a bs =1.672 in the spectrum of Q0237-0023). When systems 
suffered from severe blending such that the above mentioned 
procedure was not possible, we just include W r upper limits. 
Fe II A2383 and Fe II A2600 are the strongest of the Fe II 
transitions, and thus when Mg II is saturated they may also 
be saturated. Therefore, we measured other weaker Fe II 
transitions (Fe II A2374 and Fe II A2587), that are less likely 
to be saturated, and thus provide more leverage on an ac- 
curate measurement of the Fe II. 

In Figures [5] to 0] we characterize the Mg II absorption 
properties of our sample. Figure [5] shows the cumulative 
equivalent width distribution of Mg II A2796 for our sam- 
ple of strong Mg II absorbers. Figure shows the binned 
rest-frame Mg II A2796 equivalent width distribution. The 
distribution of eq uivalent widths is si milar to that for the 
unbiased survey of INestor et ail <|2005h , however our sample 
has a slightly larger number of moderately strong systems 
relative to very strong systems. This small difference does 
not represent a problem for the present study, since our aim 
is to compare the same types of systems at different red- 
shifts. 

Figure [4] presents the distribution of the Mg II absorp- 
tion systems in equivalent width-redshift (W r -z a bs) param- 
eter space. The top panel shows W r (Mg n) versus redshift 
for the 87 strong Mg II absorption systems in our study. 
Throughout this paper we display moderately strong Mg II 
absorption systems (0.3 < W r < 1 .0 A) as black circles and 
very strong Mg II systems (W r > 1.0 A) as red squares. The 
bottom panel shows the binned absorption redshift (z a b s ) 
distribution of all the systems in top panel (in white). The 
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Figure 1. VLT/UVES profiles of Mg n A2796, Mg n A2803, Mg I 
A2853, Fe II A2374, Fe II A2383, Fe II A2587, and Fe II A2600 
(if covered) for the strong Mg II absorption lines systems in our 
sample. The spectra are normalized and the different transitions 
are aligned in velocity space (zero point corresponds to the optical 
depth mean of the Mg II A2796 profile). The vertical dashed lines 
delimit the regions in which Mg II A2796 is detected. [See the 
electronic edition of the Journal for Figs 1.2 to 1.87.] 
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Figure 2. Cumulative Mg II A2796 equivalent width distribution 
for the systems detected in this study (shown as a solid line). 
Dashed li ne represents t he dis tribution of systems in the larger 
sample of iNestor et al.l J2005l> using the median of our sample 
((*)= 1.268). 




2 3 
W r (2796) (A) 

Figure 3. Distribution of the rest-frame equivalent width (W r ) 
of Mg II A2796 for the systems detected in this study. Most of 
the systems found lie in the lower part of the W r range. The very 
strong systems comprise ~33% (29/87) of the total sample. Only 
~ 8% (7/87) of all systems present W r > 2 A. 



87 systems are roughly uniformly distributed over the ab- 
sorption redshift range, from z a bs = 0.238 to 2.464, with 
their mean at (z) — 1.297. The redshift histogram for only 
the very strong Mg n systems (red/shaded) is superim- 
posed. The distributions in redshift for moderately strong 
((z) = 1.322) and very strong systems ((z) = 1.247) are not dif- 
ferent at a statistically significant level (the K-S tests shows 
P(K-S)=0.83 that they are drawn from the same distribu- 
tion). However, we note that there is only one very strong 
system with z > 2. 
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Figure 4. Redshift distribution of the strong Mg II A2796 in this 
study. The top panel shows rest-frame equivalent width of Mg II 
versus the absorption redshift. Black circles represent the mod- 
erately strong absorption systems (0.3 < W r < 1.0 A) and red 
squares are used for very strong systems (W r > 1.0 A). Errors in 
W r are smaller than the symbol sizes. The bottom panel shows 
a histogram of the redshift distribution. The total histogram rep- 
resents the distribution of all systems (all symbols in top panel), 
and the red area the distribution of very strong systems with 
Wr(Mg II) > 1.0 A (red squares in top panel). 



3 RESULTS 

As mentioned in !JU the Fe/Mg abundance ratio is a tracer 
of the star formation history of a galaxy. For typical physi- 
cal properties of interstellar gas, most of the Mg and Fe are 
in the form of Mg II and Fe II, respectively. Thus, although 
we detected Mg I (F e I is almost never present, however, 
sec Jones et alJ boiOf), hereafter we will focus on Mg II and 
Fe II. For some of the strong absorbers in this study, and 
more likely for the very strong absorbers, Fe II A2383 and 
Fe II A2600 are saturated, as well as Mg II (see Figure set[T|- 
Thus for these Fe II transitions, there is a slower Wr(Fe n) 
increase, with increasing Fe II column density, reducing the 
leverage of the W r (Fe n)/W / r (Mg n) ratio in determining 
evolutionary trends. We compensate for this by using weaker 
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Figure 5. Ratio of W r (Fe II) to Wr(Mg II), for the four Fe II transitions measured in this study, versus equivalent width of the blue 
member of the Mg II doublet (Wr(Mg n)). Black circles represent moderately strong Mg II absorption (0.3 < W r < 1.0 A) and red squares 
very strong Mg II absorption (W r > 1.0 A). Arrows indicate upper limits (for blends in Fe II as well as for cases with no detection), and 
error bars with no central symbol indicate cases where we estimated an upper and lower limit value for the Fe II measurement (see jj2{ . The 
number of data points is not the same among the four panels due to the non-uniform coverage of the Fe II transitions in the VLT/UVES 
spectra. As the values of WV(Mg II) increases, Wr(Fe n) values tend to increase as well, and when both profiles show signs of saturation, 
the ratio W r /W r (Fe II) approaches unity. The figure shows how Fe II A2374 and Fe II A2587 display a larger range of W r /W r (Fe II) values 
than Fe II A2383 and Fc II A2600, because the former two transitions are weaker than the latter. 



Fe II transitions, such as A2374 and A2587, which are less 
saturated, and for many systems do not show evidence of 
saturation. Since Fe n A2587 is more often covered in this 
VLT/UVES sample, we often emphasize it in drawing con- 
clusions. 

Figure [5] shows the distribution of ratios of 
W r (Fe n)/W r (Mg n) with respect to W r (Mg n) for 
the four Fe II transitions measured in this study: Fe n 
A2374, Fe II A2383, Fe II A2587, and Fe II A2600. The ratio 
of W r (Fe Il)/W r (Mg n) approaches unity as W r (Mg n) 
increases, since both the Fe II and the Mg II are saturated 
for the strongest systems. Fe II A2374 and Fe II A2587 
show the largest range of ratio values for the very strong 
absorbers since they are less saturated than the other two 



Fe II transitions. Although our number statistics are small 
for values of W r (Mg n) > 2.5A, all of our data points in 
that range have large ratios, fairly close to 1, even for Fe II 
A2374 and Fe II A2587. 

Figure |S] shows the distribution of ratios of 
W r (Fe ll)/W / r (Mg n) in absorption redshift (z a bs), again 
for the four Fe II ions measured in this study. We can see 
that these distributions for the two different absorption 
classes (moderately strong and very strong) show different 
trends. While moderately strong Mg II absorption systems 
(black circles) present a wide range of VF r (Fe ll)/W / r (Mg n) 
ratios, independent of redshift, very strong (red squares) 
Mg II absorbers are preferentially present in certain regions 
of the Wr-z a b a parameter space. For very strong Mg II 
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Figure 6. Ratio of W r {Fe ll)/lVr(Mg II) versus absorption redshift (z a (, s ) for the four Fe II transitions studied. Symbols are equivalent to 
those in Figure[5] The distribution of moderately strong Mg II absorbers seems to cover a large range of W r (Fe II) /W r (Mg II) at all rcdshifts. 
On the contrary, very strong Mg II absorbers at lower redshifts (z a f, s < 1.2) appear to cluster at high values of W r (Fe Il)/Wr(Mg II), 
which they do not do at higher redshifts. 



absorbers (red symbols) at lower redshifts (z a fcs < 1.2 in our 
sample, see below) W r {Fe n)/W / r (Mg n) is concentrated at 
higher values, and there are few absorption systems with 
small ratios of W r (Fe ll)/VVV(Mg n). As redshift increases, 
a wider range of the ratio is found. As we showed in Figure 
13] there are a sufficient number of very strong systems at 
low redshift (as compared to high redshift) that we can 
be confident that this trend is not due to small number 
statistics. 

The difference in the distribution of 
W r (Fe u)/W r (Mg n) between low and high redshift, 
for very strong absorbers, is clear from visual inspection 
of Figure [6] at low redshift, the very strong Mg n ab- 
sorbers show near-unity W r ratios while there is a much 
larger spread at high redshift, extending to smaller W r 
ratio values. We proceed to quantify the significance of 
this trend by using several statistical techniques (Table 
[2j. We investigate the probability that two subsamples 



were derived from the same population by performing 
Kolmogorov-Smirnov (K-S) and Anderson-Darling (A-D) 
tests. The A-D is a modification of the K-S test and gives 
more weight to differences in the tails of the distribution 
of the data (|Scholz fc Stephens! Il987t) . Table [2] includes 
the comparisons between the W r ratio distributions of 
moderately strong and very strong systems, and between 
high and low redshift subsamples of each. To divide each 
sample into low-redshift and high-redshift subsamples we 
compute the median of the z a bs values, excluding upper 
limits, for each Fe n transition and average them, which 
yields < z abs >= 1.20. 

The results in Table [2] confirm that the 
W r (Fe ii)/Wr(Mg n) distribution for moderately strong 
absorption systems is not consistent with that for the 
very strong absorption system population. We see that 
this difference is more significant in the low redshift bin 
(see second row of Table [2| with W, (Fe ii)/W r (Mg n) 
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approaching unity for the very strong systems (see Fig. [6]). 
Through bootstrapping, we found that in only < 0.8% of 
the realizations of the moderately strong sample at low 
redshift were there as many large W r (Fe n)/W / r (Mg n) 
values (larger than the median of W r (Fe ll)/W r (Mg n)) 
as in the very strong sample at low redshift (see Table [5] 
for more details). At higher redshifts, there is still some 
difference between the W r (Fe n)/W / r (Mg n) values for very 
strong and moderately strong samples (see third row of 
Table [2}, but it is not as pronounced. 

A comparison of the low and high redshift subsamples 
of very strong absorption systems shows that there is only 
a small probability that they share the same parent popu- 
lation (see fourth row of Table [2} . The largest probability 
of W r (Fe ll)/W / r (Mg n) of the low and high redshift sub- 
samples coming from the same parent population is for Fe II 
A2374, P(K-S)=0.04. Moderately strong absorption systems 
show no such difference (see fifth row of Table . Thus we 
conclude that there is a significant evolution in redshift in 
the W r (Fe ll)/W r (Mg n) distribution for very strong Mg II 
absorbers, but not for moderately strong absorbers. 

The left-bottom quadrant of each of the Figure [6] panels 
is almost, but not completely, devoid of very strong systems 
(red squares), thus we would like to comment on the excep- 
tions that do lie in this quadrant. Also, there are a few upper 
limits on the Fe II measurements of each transition (red ar- 
rows), due to significant contamination, and some of these 
values could, in principle, come down into this quadrant if 
they could be accurately measured. In the case of the system 
at z a b 3 — 0.8919 towards Q0300+0048 the spectrum is very 
noisy (S/N < 5; see Figure [TJ around the Fe II A2374, Fe II 
A2383, and Fe II A2587 transitions, so a point does not ap- 
pear on those panels. However, the quality of the spectrum 
at the position of Fe II A2600 is higher, and it does appear to 
yield a relatively low value of 14 7 r (2600)/W, (Mg n), which 
is apparent in the bottom-right panel of Figure A differ- 
ent situation applies in the case of the z a b s = 1.0387 system 
towards CTQ0298, another of the low redshift data points 
with a low VK r (Fe n)/W / r (Mg n) value in the bottom-right 
panel. The Fe II absorption for the transitions shortward of 
2600 A lie in the Lya forest of CTQ0298, and the resulting 
blends only allow us to place upper limits on W^Fe n) for 
those transitions. In the case of the z a bs = 0.7261 system 
toward Q0453-0423, Fe II A2374 and Fe II A2383 appear to 
be also blended with other ions, but it is likely that their 
actual W r values are close to the upper limits, which im- 
plies that their W r (Fe Il)/W r (Mg n) values are truly in the 
upper- left quadrant. We have good estimates for the un- 
blended Fe II A2587 and Fe II A2600 measurements which 
show large values of the W r (Fe ll)/W r (Mg II). Finally, all 
of the Fe II transitions associated with the z a t, s = 1.1335 sys- 
tem toward Q1621-0042 suffer from severe blending in the 
Lya forest, thus it is impossible to determine the true W r 
values. Having considered these points, we conclude that, 
with a few possible exceptions, there is a true deficit of small 
Wr.(Fe ll)/W / r (Mg n) values for very strong Mg II absorbers 
at low redshift. 



3.1 Mg II profile evolution 

Figure [6] shows there is an evolution of the 
W r (Fe ll)/14 7 r (Mg n) ratio of very strong Mg II ab- 
sorbers. We investigate the profile shape evolution of Mg II 
for the these systems in order to consider how it might 
affect the W r (Fe u)/W r (Mg n) ratio. 

Figure [7] includes all the normalized profiles of very 
strong Mg II systems in our sample, and their corresponding 
Fe II A2587, ordered by increasing z a t s . In cases where Fe II 
A2587 is not covered, another transition is used as indicated 
on the panels in the figure. We find that the profile shape 
of Mg II varies from a predominance of "boxy" profiles at 
low redshift, to the inclusion of more kinematically extended 
and less saturated Mg II profiles at high redshift, as would 
be expected for outflows/superwinds. 

Figure [5] places the very strong Mg II and Fe II A2587 
absorption profiles at their corresponding points of the 
Wr-(Fe ll)/Wr(Mg ll)-z a bs parameter space. This allows us 
to examine the relationship between profile shape and the 
W r (Fe n)/W r (Mg n) ratio. We can see from Figure [8] that 
the evolution in the Mg II profile kinematics with redshift 
seen in FigureQis related to the evolution seen in Figure|6]of 
the values of W r (Fe ll)/VK r (Mg n) of the very strong Mg II 
systems. Most of the kinematically complex, broad profiles 
correspond to small Wr(Fe ll)/W / r (Mg n) ratios at high red- 
shift. We discuss this effect in Sec i]4.1l and. in Sec i]4.2l con- 
trast it with the previous analysis of Wr(Fe ll)/W / r (Mg n) 
for weak Mg II systems. 

3.2 Kinematical properties 

To quantify the evolution we describe in section |3~T1 we also 
perform kinematic analysis on the Mg II profiles of the sys- 
tems in our sample. We saw in Figures [7] and [5] that systems 
with larger kinematically spreads tend to be more predom- 
inant at higher redshifts. We quantify this effect by calcu- 
lating the velocity spread, which is a measure of the full 
absorption line width (in units of velocity), defined as 

AW = Vlow — Vhigh, 

where vi ow and v^igh correspond to the minimum and max- 
imum velocity of the absorption, defined by a 3a absorption 
detection. Figure [5] shows the velocity spreads (Av) of the 
Mg II A2796 profiles in our sample. The distribution of Av 
for moderately strong systems (black circles) is similar at all 
redshifts, but we find differences between the very strong ab- 
sorbers (red squares) at low and high redshift. The median 
of the low and high redshift values of Av differs (255 and 
334 km s _1 , respectively) and the K-S test shows there is 
a probability of only 0.03 that both subsamples derive from 
the same population. If we omit the outliers of Av > 600 km 
s _1 , the K-S test shows the two distributions to be different 
at a even larger significance level (P(K-S) = 0.009). 

Figure [10] shows the absorbed pixel fraction (APF) dis- 
tribution of the Mg II A2796 absorption profiles, as a func- 
tion of redshift. The APF is defined as the fraction of pixels 
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Figure 7. Normalized Mg II and Fe II profiles of the very strong Mg II systems in our sample, in increasing order of z a b s , which is indicated 
at the top of each spectrum. Mg II A2796 and Fe II A2587 are plotted unless otherwise noted. There is no spectral coverage of Fe II in the 
spectrum of Q1418-064 for the absorber at z a i, a =1.4578. Velocity ranges are (-250,250) km s" 1 in the single boxes and (-500,500) km s" 1 
in the double boxes, thus the velocity scale is the same in all boxes, except for the case of Q0551-3637 at z a i, s = 1.9609, where the range 
is (-600,400) in order to include the entire profile. Dotted horizontal lines indicate flux zero. 



with detected absorption ov er the entire velocit y spread of 
the profile (Aw), as used in iMshar et~afl (120071 ). The APF 
equals one in profiles where there is absorption over the en- 
tire velocity range, and it helps to quantify the separation 
between components in the absorption profile whenever the 
values are smaller. In Figure[l0]it can be seen that there is a 
large concentration of values with continuous coverage (APF 
equal 1). The majority (70%; 61/87) of the absorbers in our 
sample show a complete absorbed coverage range, and 80% 
(70/87) show a complete or close to complete (APF > 0.95) 
coverage. We obtain slightly larger percentages if we restrict 
our computations to the very strong absorbers alone (72% 
and 86%, respectively). 



As in IMshar et all (|2007l ). we find differences between 
the low-redshift and high-redshift systems. High-redshift 
systems tend to show APF values closer to 1, which indi- 
cates that the absorption cover the full velocity range, while 



at low redshift it is more likely to find APF values lower than 
1. This is due to a tendency for the high-redshift systems to 
have weak components connecting the stronger absorbing 
regions. 

Because in the case of very strong absorbers, so many 
APF values are close to unity, this measurement does not 
provide much leverage in determining kinematic evolution. 
To further quantify the shape of the profiles, and in par- 
ticular how "boxy" or spread out they are, we define their 
dimensionless average depth as the fraction of normalized 
flux that is absorbed, averaged over the full range of Av. 
In the case of a completely rectangular absorber, the depth 
would be equal to one (since our W r is defined from normal- 
ized spectra). This averaged dep th is equivalent to the D- 
in dex as defined in Ellison! l|2006l ). and posteriorly reviewed 
in lEUison et all (2009), except that we include the whole 
velocity range (An) in the average. We present this depth 
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Figure 8. Mg II A2796 and Fe II A2587 profiles of very strong Mg II systems (W r > 1.0 A) in the parameter space of 
log(W r (Fe n)/W / r (Mg II)) vs. z, which was previously displayed in Fig[6] (symbols are equivalent). For ease of location, we have grouped 
contiguous systems in boxes and circles. The figure shows that while "boxy" profiles appear at both low and high redshift, they are more 
predominant at low redshift. On the other hand, systems with "outflow/superwind" type profiles (less saturated and showing larger velocity 
spreads) are only present at high redshift, and correspond to the systems with smaller values of Wr(Fe ll)/Wr(Mg n). 
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Figure 9. Velocity spread versus absorption redshift of the 87 
systems analyzed. Symbols are equivalent to those in Figure \5\ 
Strong Mg II absorption profiles (black circles) show a homoge- 
neous distribution with redshift, but very strong Mg II absorp- 
tion profiles (red squares) show a significant difference in their 
median value of velocity spread between low (z a {, s < 1.2) and 
high (z a i, s > 1.2) redshift. 



measurement versus some other properties of our sample in 
Figures 1111 1121 and 1131 Figure [11] and Figure [12] show pro- 
file depths versus W r (Mg n) and Av, respectively. There 
is an envelope at the lower end of the depth distribution 
which depends on the W r (Mg n) value. This is regulated 
by the typical maximum Av value for the profiles. In other 
words, we do not find shallow profiles with large W r (Mg n) 
values (see Figures [7] and [8]l , which would result from very 
large Av values and would lead to a departure from the 
W r (Mg ll)-Av correlation. Only the several profiles with 
the largest W-^Mg n) have Av > 500-600 km s" 1 . Thus, a 
given equivalent width can only result from a depth larger 
than some minimum value. Completely "boxy" profiles are 
also rare, since weaker subsystems are common; because of 
this and the gradual recovery of profiles to meet the contin- 
uum, there is also a natural high limit for the depth value. 
Figure [12] shows that moderately strong and very strong 
Mg II absorbers occupy a different region of the depth-Au 
parameter space. This is also partly due to the correspon- 
dence between larger Wr and larger Av. Once saturation is 
present the only way to increase the W r is through increas- 
ing the Av of the profile. 

Figure [13] shows the distribution of depth with absorp- 
tion redshift. There is a clear evolution for the very strong 
Mg II absorption profiles, which is not present for the mod- 
erately strong ones. An Anderson-Darling test results in a 
probability of only P = 0.03 that the depths of the low 
redshift subsample of very strong absorbers are drawn from 
the same population as those of the high redshift subsam- 
ple. The distribution of depths of the moderately strong 
(black circles) and very strong absorbers (red squares) at 
low redshift is significantly more different (P(K-S)= 4.e-7) 
than the distribution of the two subsamples at high redshift 



Figure 10. Absorbed pixel fraction versus absorption redshift of 
the systems in our sample. Values of 1 indicate systems where the 
whole velocity range is absorbed. Symbols are equivalent to those 
in Figure [5] 



(P(K~S)— 0.02). Because the depth is influenced both by 
W r (Mg n) and Av, evolution of W r (Mg n) would play a 
role into depth evolution. However, in Figure [4] we already 
showed that there is no significant evolution of W r (Mg n) 
with redshift in our sample, thus changes in Av must domi- 
nate. Another way to consider the trend of depth evolution 
in Figure [13] is to note that the cluster of very strong ab- 
sorbers with depth values ~ 0.4 in Figure [121 are all in our 
high redshift sample. The presence of small depths at high 
redshift, but not at low redshift, is reminiscent of the evo- 
lution of W r (Fe ll)/W,-(Mg n) shown in Figures [6] and we 
will consider this connection in § [4] 



4 SUMMARY AND DISCUSSION 

In a large sample (81) of high resolution VLT/UVES spec- 
tra, we have found 58 moderately strong (0.3 < W r < 1 A) 
and 29 very strong (W r > 1 A) Mg II absorption systems, 
and searched for accompanying Fe II and Mg I absorption. 
We have investigated the profile shape evolution of Mg II 
as well as the evolution of the Wr-(Fe ll)/W r (Mg n) ratio. 
Here we summarize the main results of this study. 

• For moderately strong Mg II absorbers, the ratio 
W r (Fe u)/W r (Mg n), does not evolve significantly from 
z ~ 2.5 to z ~ 0.2. 

• On the contrary, very strong Mg II absorbers show evo- 
lution in W^Fe n)/W r (Mg n), with a deficit of small values 
at low redshift. 

• Moderately strong Mg II absorber profiles do not evolve 
significantly from z ~ 2.5 to 2 ~ 0.2, as quantified by profile 
velocity spread and depth. 

• However, very strong Mg II absorbers at low redshift 
all have relatively high profile depth values, while those at 
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Figure 11. Profiles' depth distribution with restframe equivalent 
width of Mg II A2796. Depth is defined as the fraction of flux ab- 
sorbed, averaged over the velocity width of a profile. Symbols are 
equivalent to those in Figure [5] Moderately strong Mg II absorp- 
tion profiles (black circles) show a larger range of depths than 
very strong Mg II systems (red squares), which show a minimum 
depth of ~0.3. This is not surprising because there is a correlation 
between W r (Mg ll) and Av. 



Figure 13. Depth versus absorption redshift of the profiles in our 
sample. Symbols are equivalent to those in Figure [5] We can see 
that while the moderately strong Mg II absorption profiles (black 
circles) show a similar distribution at low and high redshift, very 
strong ones (red squares) show evolution with a larger range of 
depths at high redshifts than at low redshifts. 
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Figure 12. Depth versus Av of the profiles in our sample. Sym- 
bols are equivalent to those in Figure[5] Moderately strong (black 
circles) and very strong (red squares) Mg II absorbers are confined 
in different regions of the depth-Au parameter space. 



high redshift have a wider range of values (more similar to 
the range for moderately strong systems). 



4.1 Possible Causes for the Evolution of Very 
Strong Mg n Absorbers 

It is not surprising to find an evolution in the p opulation 
of very strong Mg II systems. iNestor et alj |2005l ) already 
showed an excess of these systems at high redshift, rela- 
tive to absorbers with smaller equivalent widths. Among the 



very strong Mg II absorber population, the deviations from 
cosmological evolu tion are larger as W r {Mg ll) increases 
ijNestor et al.ll2005t ). 

Given the observed evolution in the 
W r (Fe Il)/VFr(Mg n) ratio for the class of very strong 
Mg ll absorbers, shown in Figure [8] we consider several 
plausible reasons. These include: 1) a changing ionization 
rate due to evolution of the extragalactic background; 2) 
a change in the level of a-enhancement due to changing 
contributions of Type II and Type la supernovae to the 
absorbing gas; 3) a change in the kinematics of the profiles, 
including different levels of saturation and velocity spreads 
which affect the equivalent widths of Fe II and Mg II 
differently. Here we will consider the expected effect of each 
of these factors. 

The parameter space of very strong Mg II absorbers 
shows large W r (Fe n)/W / r (Mg ll) at both low and high red- 
shift, but small W r (Fe Il)/W r (Mg n) only at high redshift 
(see Figs. [6] and [Sj . The profiles of very strong Mg II ab- 
sorbers (those with W r >1 A) are almost always saturated 
in Mg ll, and sometimes saturated even in the weakest ac- 
cessible Fe ll transition as well (see Fig. 0. We should first 
consider this when interpreting the observed absence of low 
W r (Fe ll)/W r (Mg ll) values in this population at low red- 
shift. As shown in 13.11 Figure [8] indicates that there is an 
evolution in the Mg II profile kinematics that is relevant for 
our interpretation of the W r (Fe ll)/W / r (Mg ll) evolution for 
very strong Mg II absorbers. At low redshift, many of the 
systems show "boxy" profiles and they tend to have fewer, if 
any, satellite clouds around their compact, "boxy" portions. 
Therefore, they do not have a significant range of velocity 
that is unaffected by saturation in Mg II. A typical example 
of these "boxy" profiles is the z — 0.2378 system toward 
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Q0952+0179, shown in Figure [Q (also see Fig.[7]and Fig.©. 
Because the weaker Fe n A2374 and Fe II A2587 transitions 
are not strongly saturated for many components, the Fe n 
equivalent widths are less affected by saturation effects. 

Also, it is important to note that we are only com- 
paring the equivalent widths of the Mg n and Fe n tran- 
sitions, and not their column densities. In the case of a sat- 
urated profile, of course it is not possible to measure the col- 
umn density, and we should interpret the equivalent width 
as a lower limit on how much material is present. For the 
W r (Fe ll)/W r (Mg n) ratio, in cases where the Mg n is sat- 
urated and the Fe n is not, our values are therefore up- 
per limits. This alone could lead to the absence of very 
strong Mg n absorbers with small W r (Fe Il)/W r (Mg n) 
values. The question then is why very strong Mg n ab- 
sorbers at high redshift do often have small values of 
W r (Fe ll)/W r (Mg n), while those at low z do not. 

The answer can be seen by examining the Mg II pro- 
file shapes of the high redshift systems in the lower right 
quadrants of Figures [6] and [8] These profiles show compo- 
nents more spread out in velocity space, with several subsys- 
tems within the absorption profile. Because of this kinematic 
spread, even if the column densities and numbers of compo- 
nents are the same as for the low redshift systems, we will see 
less saturation. This is evident in the existence of some low 
values of profile depth (see Figure [13)) among the high red- 
shift very strong Mg n absorbers. Several "boxy" saturated 
Mg n profiles are apparent among the high z systems as well, 
and these have large values of VF r (Fe ii)/W r (Mg n), as did 
all of the systems at lower redshift. If the Mg II profile suffers 
from significant saturation, the W / r (Fe n)/W / r (Mg n) value 
is an overestimate and it is impossible for it to be small. In 
order that a very strong Mg II absorber not have significant 
saturation it must be kinematically spread, and such kine- 
matics are apparently more common at high redshift than 

at low redshift. 

This is consistent with the conclusions of iMshar et al.l 

(120071 1 on the evolving kinematics of a smaller sample of 
very strong Mg II absorbers. That study found that complex 
systems with many components spread in velocity are more 
common at high redshift than at low redshift. F or example, 
profiles typically associa ted with superwinds (|Bond et al.l 
l200ll ; lEllison et alJlioblh show large velocity spreads. The 
ex cess of very strong M g 11 absorbers at high redshift, found 
in lNestor et al.l |2005h . may be related to our findings if the 
excess systems correspond to outflows/superwinds, and are 
thus evolving together with the star formation rate in galax- 
ies. 

Returning to the three effects that could affect evolution 
of the W^Fe ii)/W r (Mg 11) ratio in very strong Mg 11 sys- 
tems, we have found that kinematic evolution of this popu- 
lation could explain the observed trend. This does not neces- 
sarily mean that ionization and abundance pattern changes 
are not taking place as well, just that we are not able to 
get an indication of the ionization parameter or abundance 
pattern for the saturated systems. In other words, some of 
the systems with "boxy" profiles in Mg 11 both at high and 



low redshift, could have higher/lower ionization parameters 
or could be a-enhanced/depleted relative to the others, but 
would still not have small W r (Fe ll)/W / r (Mg 11) ratios rela- 
tive to others because of saturation of the Mg II. At higher 
redshift, however, it seems likely that systems with large 
kinematic spread and small values of W r (Fe ll)/W / r (Mg 11) 
(bottom right quadrant of Figure [8} are different than the 
"boxy" profiles systems (top right quadrant of Figure [8} 
that also appeared in this redshift range. The strong Mg 11 
absorbers, with 0.3 A < W r (Mg 11) < I A, also have the 
full range of Wr-(Fe n)/W r (Mg 11), but not all of them with 
large values are strongly saturated in Mg II. 

It is natural to expect an absorber population to be 
more ionized at higher redshift, because of the evolution of 
the extragalactic background radiation. This is due both to 
the changes in the ionizing extragalactic background radia- 
tion and due to the more intense local radiation field in star- 
bursts, both of which are a consequence of the higher global 
star formation rate at higher redshift. Even without a con- 
tribution from local star formation, the photon number den- 
sity, logn 7 changes from —4.7 cm -3 at z = 2 to —5.7 cm~ 3 
at z — 0.3, which will lead to a significant difference in the 
ionization parameter affecting our Mg 11 absorbing clouds. 
In particular, JV(Fe ll)/iV(Mg 11) is relatively constant for 
ionization parameter log U < —4.5 but decreases rapidly 
from log U — —4 up to higher values. At z — 2, a cloud 
with Tin = 0.1 cm -3 would have log U = —3.7, while it 
would have log U = —4.7 at z = 0. If instead the Mg 11 
cloud has lower densities, for example uh = 0.01 cm -3 , it 
would result in log U = —2.7 at z — 2, and log U = —3.7 at 
z = 0. While the difference of iV(Fe n)/7V(Mg 11) between 
low and high redshift systems would be barely noticeable 
in the case of ng =0.1 cm' 3 , clouds of densities of uh = 
0.01 cm -3 would have significantly lower N(Fe n)/AT(Mg 11) 
at high redshift than at low redshift. Thus, the systems 
with small N(Fe ll)/iV(Mg 11) values would result from sys- 
tematically lower density clouds. A larger range of densi- 
ties (and thus log U) would lead to ionization of the Fe 11 
and would give rise to small N(Fe ll)/iV(Mg 11) at high red- 
shift, and thus we would expect more systems with small 
N(Fe ll)/7V(Mg 11) at high redshift. Although this clearly 
matches the observed behavior for very strong Mg 11 ab- 
sorbers, we cannot measure this effect since it is disguised 
by the saturated Mg 11 profiles in many of the systems. 

Similarly, a-enhanced systems are expected to be more 
common at higher redshift, where the star formation rate 
is higher. It takes time for a stellar population to produce 
Type la supernovae, leading to about a 1 billion year de- 
lay before the iron abundance would be elevated to the so- 
lar value, re lative to magnesi u m. A s will be discussed fur- 
ther in 14.21 iNaravanan et all (|2007l ) noted that this delay 
could be a significant contributing factor to the absence 
of weak Mg 11 absorbers with large W r (Fe n)/W r (Mg 11). 
The larger kinematic spreads of some of the high-redshift 
very strong Mg 11 absorbers may be related to enhanced 
star formation/starburst and wind activity at high redshift. 
Thus, these systems would be expected to be the a-enhanced 
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systems with small W r (Fe n)/W / r (Mg n) values. It is plau- 
sible that there would be fewer such a-enhanced systems 
at low redshift. Also, a stronger depletion of Fe relative 
to Mg would result in a smaller W r (Fe ll)/W r (Mg 11). It 
is known that Mg returns to the gas phase faster than Fe 
l|Fitzpatrickl 1 19961 ) . Indeed, mild and high Fe depletion has 
been previously observed in some DLA and sub-DLA sys- 
tems, both at low and high redshift (e.g .. Ledoux et alj|2002l; 
Meiring et~all |2008| ; iNoterdaeme et al.|[2oTol ; iMeiring et al l 



201 if ). Like ionization effects, both a-enhancement and de- 



pletion of Fe contribute in the right direction to lead to a 
relative absence of very strong Mg absorption in the lower 
right quadrant in the W / r (Fe n)/W / r (Mg n) vs. z diagram. 
However, we cannot tell if the lower left quadrant is vacant 
because of the lack of a-enhancement/greater depletion of 
Fe at low redshift, or just because low redshift systems pref- 
erentially have kinematics that lead to saturation of Mg II 
and Fe II. 

We conclude that, in the case of very strong Mg II ab- 
sorbers, low values of W r (Fe n)/W / r (Mg n) at high redshift 
could be due to a larger ionizing radiation field or to a- 
enhancement of a larger fraction of systems, but the effects 
cannot be quantified in these mostly saturated systems. We 
also cannot be certain of the cause of the absence of low 
W r (Fe ll)/W / r (Mg n) values at low redshift, because all the 
low redshift systems in our sample have significant satura- 
tion of Mg II and Fe II. The significant evolution that we 
have detected, therefore, is the evolution of the kinematics 
of the Mg II profiles, from a mixture of saturated, flat bot- 
tom profiles and kinematically spread components at high 
redshift, to only predominantly saturated, flat bottom pro- 
files at low redshift. 

This is quite consistent with several recent studies of 
the nature of DLAs and sub-DLAs, which are coincident 
with the population of very strong Mg II absorbers. Since 
very strong Mg II absorbers are selected by an equivalent 
width criterion, they are not necessarily all among the high- 
est column density absorbers, and thus are not all DLAs. 
The equivalent width of Mg II can be very large either be- 
cause the total column density of material is very large (thus 
it would likely be a DLA), or because the kinematic spread 
of components is large enough that the equivalent widths 
add up to a large total (in which case the system might be 
a DLA, but may also be a sub-DLA or even a Lyman limit 
system) . 

Indeed, the populations of DLAs and sub-DLA s have 
important intrinsic differences. iKulkarni et al.l (|2010h show 
that DLAs have systematically lower metallicities than 
sub-DLAs (based on Zn measurements), and furthermore 
that they have a sm aller velocity spread than sub-DLAs. 
IKulkarni et al.l (|2010T ) suggests that the differences in metal- 
licities between DLAs and sub-DLAs could be due to dif- 
ferent star formation histories, where the galaxies observed 
as sub-DLAs underg o m ore rapid star form ation. Both 
IKulkarni et al.l (|2010T) and lNestor et all l|201ll ) suggest that 
DLAs are often dense regions in ordinary galaxies, while 
sub-DLAs are of higher metallicity and extended kinemat- 



ics such as one would expect from a superwind. All these 
results would be consistent with the boxy profiles that 
we see in Mg II tending to be DLAs, while the kinemat- 
ically spread systems that are more common at higher 
redshift, would tend to be sub-DLAs, where the star for- 
mation is more likely to happen. Then, our finding that 
small W^Fe ll)/VK r (Mg n) values for very strong Mg II ab- 
sorbers are only present at high redshifts could be due to 
both a-enhancement and kinematics, expected from out- 
flows /superwinds . 

Another related factor that would condition the shape 
of the Mg II profiles may be the association of the subset of 
very strong absorbers with galaxy over-density regions such 
as group s or clusters where interactions and mergers are 
common. iNestor et all (|2007l ) found that field imaging of a 
subset of the strongest Mg II absorbers (W r (Mg n) > 2.7 A) 
at low redshift (0.42 < z < 0.84) indicates that pairs of 
galaxies, often with distorted morphologies, and starburst 
related phenomena, are coincident with these Mg II sys- 
tems. Although, this redshift range is lower compared to 
the high-2 very strong systems in our sample, interactions 
would be only more frequent at high redshift. If these very 
strong Mg II absorption systems are selecting group envi- 
ronments, star formation is also going to be higher com- 
pared to the average for t hat redshift (e.g. Kennicutt et all 
1 19871 ; and more recently IWong et al.l l201ll ). which would 
push this subset of absorbing gas in the direction towards 
lower W r (Fe n)/W r (Mg n) ratios. We find that our results 
are consistent with this scenario. In our sample, for one 
of the stron g est M g n systems (Q0002-422, z = 0.837), 
lYannv et al. ( 1990) have identified several emission line 
sources within 100 kpc, which indic ates that there is prob- 
ably a group environment (also see lYork et alj Il99ll ) . Sim- 
ilarly, the very strong Mg II system in the line of sight of 
Q0453-423 a t z = 0.726 also h as multiple galaxies associ- 
ated with it l|Yannv et al.lll990l ') ■ 



4.2 Comparison to the Evolution of Weak Mg II 
Absorbers 



The evolution of 
W r (Fe u)/W r (M g n) 



weak 
ratio 



systems and their 
was present ed 



m 



INaravanan et al.1 ((2007) and iNaravanan et all (2008). 
Figure Q3] shows the evolution of the W r ratio for very 
strong Mg II absorption systems (red squares - top) and 
the weak Mg H s ystem s (green triangles - bottom) from 
INaravanan et aT, | 120071 ). For very strong absorption (top) 
IFr(Fe n) is measured in the Fe II A2587 transition because 
it is weak enough not to be saturated, but is the most 
often observed Fe II transition in our sample. Weak Mg II 
absorption is accompanied by even weaker Fe II troughs, 
and since Fe II A2383 is the strongest transition, it is the 
mostly likely to be detected and the one displayed in Figure 
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Figure 14. Ratio of W r (Fe ll)/W r (Mg II) versus absorption red- 
shift (z a (, s ) for the very strong Mg II systems (top) and weak 
Mg II systems (bottom). Symbols are equivalent to those in Fig- 
ure [5] Two opposite trends arc present. There is an absence of 
very strong absorbers with low ratios at low redshift (top) and 
an absence of weak absorbers with large ratios at high redshift 
(bottom). In the case of the very strong absorbers, we chose 
Fc II A2587, because it was less often saturated than the Fe II 
A2383 and A2600 transitions, and was covered more often than 
the weaker Fe II A2374 transition. For weak Mg II absorbers, Fe II 
A2383 was chosen because, as the strongest Fe II transition, it is 
mostly likely to be detected. 



these weak profiles is fairly rare, the cause of the evolution 
could be 1) changing ionization or 2) abundance pattern 
changes in the population of gas clouds that produce these 
absorbers. In the first case, because of the known evolution 
of the extragalactic background radiation (higher at higher 
redshifts) and since Fe 11 is more readily ionized than Mg II, 
we would need higher density absorbers to populate the large 
W r (Fe ii)/Wr(Mg 11) high redshift quadrant, than we would 
at low redshift. In the second case, a delay between the Type 
II and Type la SNe, i.e., an a-enhancement, might help 
explain the absence of large W r (Fe n)/W r (Mg 11) values at 
high redshifts. There may not yet have been time for Type 
la supernovae to be active in weak Mg 11 absorbing clouds at 
high redshift, which themselves ten d to be less common at 
z > 1.4 than at lower redshifts (see | Naravanan et al. I l2008l ; 
iNaravanan et al.ll2007l ; iLvnch et af]|2006h . 

In fact, the evolution of the very strong Mg 11 ab- 
sorbers and the weak Mg 11 absorbers, though in the op- 
posite sense, may be related. Weak Mg 11 absorbers with 
iV(Fe 11) ~ iV(Mg 11) require an abundance pattern influ- 
enced by Type la supernovae, and their relative absence at 
z > 1.4 implies that they are related to a population of gas 
clouds that is not comm on at high redshifts. As described in 
INaravanan et all (|200ST) . this could either be because they 
are related to dwarf galaxies, in which star formation peaks 
later, or to a rarity of weak Mg 11 clouds because halos are 
more crowded so they are more commonly combined to form 
strong Mg 11 absorbers. At least some weak Mg 11 absorbers 
are produced by the same gas clouds that produce the "satel- 
lite clouds" of strong Mg 11 absorbers, simply viewed from 
a different vantage point. These are also likely to be analo- 
gous to Milky Way high velocity clouds, at various redshifts 
l|Richter etld1l2009l) . 

What we have seen here is that some fraction of very 
strong Mg 11 absorbers, those that are related to sub-DLAs 
and to star formation activity, have components spread in 
velocity space. Their W(Fe n)/W(Mg 11) values are consis- 
tent with a-enhancement, just as are all of the weak Mg 11 
absorbers in the same redshift range. Since satellites around 
very strong Mg 11 absorbers are analogs to isolated weak 
Mg 11 absorbers, and since these are common in high red- 
shift sub-DLAs, the picture of a-enhancement in both pop- 
ulations is consistent. 



I14FI Weak absorbers reflect a different trend than very 
strong absorbers: we find that large W,-(Fe n)/W r (Mg 11) 
values are uncommon at higher redshifts, i.e., weak Mg 11 
systems show a deficit as well, but in the opposite quadrant 
of t he VK r (Fe ll)/W r (Mg 11) - z parameter space. 

INaravanan et al.l (|2007l ) considered the possible causes 
of the evolution of weak Mg 11 absorbers. Since saturation of 

3 Table 1 in Narayanan et al. 2008 lists the weak absorbers 
W r (Mg 11) and W r (Fe II2383) measurements and limits. 



4.3 Conclusions and Future Work 

The main result of this study is the evolution of the popu- 
lation of very strong Mg 11 absorbers (Wr(Mg 11) > 1.0 A) 
over the redshift range, 0.2 < z < 2.5. There is a signif- 
icant absence of small W r (Fe n)/W r (Mg 11) values at low 
redshift, as compared to the full range apparent at high red- 
shift. A major contributing factor to this change is cleary 
the kinematic evolution of the profiles of the very strong 
absorbers. At low redshift, Mg 11 presents simple, boxy, sat- 
urated profiles, while at high redshift there are also Mg 11 
profiles with multiple components, some unsaturated, that 
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are spread over a larger range of velocity. Thus at high red- 
shift, it is more likely that some of the Fe II profiles will also 
be unsaturated, so that small W r (Fe ll)/WV(Mg n) values 
can arise. 

Other factors could also contribute to the absence of 
small W r (Fe Il)/W r (Mg n) values at low redshift. Most no- 
tably a small ratio could be indicative of an a-enhanced 
abundance pattern, which would be characteristic of a young 
stellar population, such as one contributing to active galac- 
tic outflows. The absence of small ratios of Fe II to Mg II at 
low redshift, is thus consistent with the idea that outflows 
are more active at high redshift. 

The key to testing the hypothesis that abundance pat- 
tern is evolving, in addition to the known evolution in the 
kinematics, is to determine if the very strong Mg n ab- 
sorbers at z > 1.2 that have small W(Fe n)/W(Mg il) are 
truly a-enhanced. This requires a separation of saturation 
effects from the effect of changing abundance patterns. Tran- 
sitions with lower oscillator strengths should be studied to 
facilitate this separation. There are numerous Fe II tran- 
sitions with various oscillator strengths in the rest frame 
A = 1050-2600 A regime. Since other Mg II transitions are 
too weak to be detected, even in these very strong systems, 
we will need to consider other a-elements for comparison, 
e.g. Si II, which has several promising transitions in the rest 
frame far-UV, an ionization potential similar to Mg II, and 
similar d ust abundances at least in the Gala ctic disk and 
the ISM (|Cartledge et al.ll2006l ; IWhittetll2010h . High resolu- 
tion quasar spectra with wide wavelength coverage (optical 
and UV) are thus essential to assessing evolving chemical 
abundance processes over cosmic time. 
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Table 1: Rest-frame Equivalent Widths for Target Transitions 
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1 


.637 





,551± 


0.002 


<0.018 


< 0.19 




0400+0 0009 

U . U4:UU_LU. UUUc? 


0.125± 


0.001 


0.049T 





.001 




1 


.657 


0. 


,683± 


0.001 


01 01 +0 0006 

U .U 1U -L _I_U. UUUU 


0.107± 


0.001 




< 0.20 




0.029± 





.001 




1 


.672 


1. 


,283± 


0.001 


956+ 001 

U.ZUU_I_ U.UU-L 


0.545± 


0.001 


354 - 496 


0.551± 


0.001 


< 0.76 






C1TO0998 


1 


.039 


1. 


,554± 


0.006 


< 0.90 


< 0.82 




< 0.37 


0.736± 


0.003 


0.176± 


0. 


.008 


O0300+0048 

v^UOUU i UU^LO 





.892 


1 


,05± 0.01 


0.18 


0.30 




0.07 


0.50± 0.01 


0.330± 





.01 


O0398-0979 





.788 





,53± 0.01 


01 6+ 003 

U.U-LU—L U.UUO 


0.21± 0.01 


057+ 009 

U.UO ( _l_ U.UUc/ 


0.109± 


0.007 










1 


.123 


1 


,42± 0.01 


409+ 008 

U.4Ui/_L U.UUO 


0.912± 


0.008 


5QQ+ 007 

U.Uf3-L U.UU t 


0.946± 


0.008 










1 


.299 





,500± 


0.008 


045+ 006 

U.U^U-L U.UUU 


0.098± 


0.007 


098+ 005 

U.UZO—L U.UUO 


0.048± 


0.005 


0.025± 





.004 




1 


.307 





,584± 


0.009 


< 0.043 


0.078± 


0.006 


093+ 005 

U.UZO—L U.UUO 


0.045± 


0.006 








00399-0385 





.763 


0. 


,615± 


0.005 


091 + 009 

U.UZjJ l_ U.UUZi 


0.156± 


0.003 


1 03+ 005 

U. -LUO— 1— U.UUO 


0.166± 


0.005 


0.055± 


0. 


.003 




1 


.438 





,385± 


0.004 








090 - 040 

U.UZU U.U4:U 


0.056± 


0.003 


0.039± 





.002 


00499-4901 





.554 





335± 


0.003 


090+ 003 

U.UZU—L U.UUO 


0.080± 


0.004 


01 5+ 009 

U.U-LO—L U.UUZi 


0.085± 


0.004 


<0.0039 






1 


.119 





,317± 


0.002 


03Q+ 009 

U.UOf— 1— u.uuz 


0.190± 


0.003 


0QQ+ 001 

U.Uc/C/— 1— U.UU-L 






0.059± 





.002 




1 


.355 





,791± 


0.002 








51 3+ 002 


0.644± 


0.002 


0.119± 


0. 


.002 


O0453-0423 





.726 


1. 


366± 


0.002 


< 0.51 


< 1.14 




710± 001 


1.052± 


0.001 


0.455± 





.002 







.908 





,853± 


0.001 


0.119± 0.002 


0.373± 


0.002 


0.186± 0.003 


0.374± 


0.002 


0.119± 


0. 


.001 




1 


.150 


4 


,461± 


0.002 


1.778± 0.002 


3.492± 


0.001 


2.693T 0.002 


3.658± 


0.001 


1.524T 





.002 




1 


.630 





,315± 


0.001 


0.0041±0.0008 


0.015± 


0.001 


0.0057±0.0008 


0.0052±0.0005 










2 


.305 





,464± 


0.001 


0.028 - 0.044 


0.154± 


0.001 








< 0.45 






Q0549-0213 





.440 





,346± 


0.007 








0.24± 0.02 


< 0.22 




0.034± 





.008 


Q0551-3637 


1 


.225 





,328± 


0.004 


0.006± 0.002 


0.079± 


0.002 








0.039± 





.003 




1 


.961 


5. 


,15± 0.03 


1.732± 0.009 


3.375± 


0.008 


2.821 - 3.659 


< 3.6 




1.090± 





.01 


Q0926-0201 


1 


.106 





,375± 0.004 


<0.0026 


0.045± 


0.002 


0.007± 0.001 


0.038± 0.002 


0.019± 





.002 


Q0940-1050 


1 


.789 


1 


,129± 0.002 


0.040± 0.001 


0.284± 


0.001 


0.095 - 0.111 


0.280 - 


0.300 








Q0952+0179 





.238 


1 


107± 0.006 








0.65± 0.01 


0.84± 0.01 


0.361T 





.008 


Q 1122- 1648 





.682 


1. 


833± 0.002 


0.399 - 0.441 


1.2786±0.0008 


0.832± 0.001 


1.3125±0.0008 


0.140± 





.001 



Continued on Next Page. . . 
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Table 1 - Continued 



QSO 


Zabs 


W T (Mg 


II) (A) 


W r (2374) (A) 


W r (2383) (A) 


W r (2587) (A) 


W r (2600) (A) 


W r (2853) (A) 


Q1127-0145 





.313 


1. 


771± 0.005 


0.89± 0.02 


1.18=b 0.01 


1.204± 0.005 


1.331± 0.005 


1.009± 0.005 


Q1151+0068 


l 


.774 


0. 


716± 0.005 


0.411± 0.005 


0.546± 0.004 


0.518± 0.004 


0.566± 0.004 






l 


.819 


1. 


07± 0.02 


0.049± 0.005 


0.255± 0.005 


0.105 - 0.235 


0.192 - 0.268 


0.138± 0.009 


Q1157+0014 


l 


.944 


1. 


539± 0.006 


0.882± 0.005 


1.157± 0.005 


< 1.1 


< 1.3 


0.716± 0.007 


Q1202-0725 


l 


.755 


1. 


899± 0.007 






0.924± 0.006 


1.279± 0.006 


0.180 


Q1229-0021 





.395 


2 


064± 0.006 






1.16± 0.01 


1.43± 0.01 


0.641± 0.007 


Q1246-0217 





.727 





453± 0.003 


0.015± 0.003 


0.142± 0.005 


0.060± 0.006 


< 0.34 


0.061± 0.003 


Q1331+0170 


1 


.777 


1. 


259± 0.003 






0.673± 0.002 


0.823± 0.001 


0.334± 0.003 




1 


.786 


0. 


878± 0.003 






0.311 - 0.369 


0.505± 0.002 


0.137± 0.003 


Q1337+0113 


1 


.637 





,617± 0.005 


0.054± 0.003 


< 0.30 


0.156± 0.004 


0.326± 0.003 




Q1341-1020 





.873 





,528± 0.006 


0.029± 0.005 


0.208± 0.006 


0.092± 0.005 


0.210± 0.004 


0.078± 0.006 




1 


.277 


1 


,454± 0.008 


0.272± 0.007 


0.802± 0.005 


0.496± 0.007 


0.812± 0.005 


0.303± 0.007 


Q1418-0064 


1. 


.458 


2 


,19± 0.01 










0.200± 0.01 




1 


.982 





,453± 0.008 


<0.0058 




<0.010 


0.011± 0.002 


0.013± 0.005 




2 


.257 


0. 


,62± 0.01 


<0.0060 


<0.0086 


<0.011 


< 0.046 


< 0.040 


Q1444+0014 





.660 





,341± 


0.005 


< 0.023 


0.096± 0.003 


0.022± 0.002 


0.086± 0.004 


0.044± 0.002 




1 


.159 





,722± 


0.004 


0.010± 0.002 


0.139± 0.003 


0.040± 0.003 




0.016± 0.002 


Q1621-0042 


1 


.133 


3. 


,168± 


0.004 


< 1.1 


< 2.3 


< 1.2 


< 3.0 


0.440± 0.005 


Q1629+0120 





.900 


1. 


,022± 


0.007 






0.383± 0.007 


0.70± 0.01 


0.205± 0.007 


Q2000-0330 


2 


.033 





,811± 


0.002 


< 0.26 


0.493± 0.001 


< 0.35 


0.490± 0.002 






2 


.292 





,327± 


0.003 


0.0121±0.0006 


< 0.32 


0.071± 0.003 




< 0.045 


Q2044-0168 


1 


.329 





,573± 


0.005 










0.032± 0.004 


Q2126-0158 


2 


.022 





,750± 


0.003 


< 0.071 


< 0.35 


0.106± 0.001 


0.282± 0.002 




Q2206-0199 





.752 





,885± 


0.004 


< 0.032 


< 0.20 






0.232± 0.005 


Q2217-2818 





.942 





,571± 


0.001 


0.0059±0.0005 


0.098± 0.001 


0.0307±0.0006 


0.0867±0.0007 


0.0288±0.0008 




1 


.628 


0. 


,656± 


0.001 


0.0040±0.0003 


0.0573±0.0005 


0.0239±0.0007 


0.0406±0.0007 


0.0348±0.0007 




1 


.692 


1 


,694± 


0.001 


0.0452±0.0007 


0.3381±0.0007 


0.1159±0.0009 


0.346± 0.001 


< 0.21 


Q2225-2258 


1 


.412 





,418± 


0.003 


<0.0021 


0.076± 0.002 


0.013± 0.001 


0.052± 0.002 






1 


.639 





,347± 


0.002 


0.004± 0.001 


0.044± 0.001 




0.039± 0.002 


0.021± 0.002 


Q2314-0409 


1. 


.045 


0. 


,349± 


0.005 


<0.0064 


0.033± 0.005 


<0.0064 


0.032± 0.005 


<0.0073 


3c336 





.318 





,48± 0.01 










<0.014 







.472 


0. 


,809± 0.009 






0.11± 0.01 


0.26± 0.01 


0.035± 0.006 







.656 


1 


,428± 0.009 


0.48± 0.01 


1.05± 0.01 


0.80± 0.01 


1.11± 0.01 


0.213± 0.009 







.797 





,46± 0.01 


<0.0089 


0.029 - 0.091 


<0.0083 


0.028± 0.003 


< 0.0078 







.891 


1 


,520± 0.006 


0.40± 0.01 




0.744± 0.007 


1.114± 0.007 


0.280± 0.006 



This table includes the measurements of the 87 strong Mg II systems detected in our sample of 81 quasar spectra. Col (1): name 
of the quasar. Col (2): redshift of the absorber. Col (3)-Col(8): measured rest-frame equivalent widths of Mg II A2796, Fe n A2374, 
Fe n A2383, Fe II A2587, Fe II A2600, and Mg I A2853, respectively. If the transition was covered, but not detected, 3cr upper limits 
are included after a "<" sign, while"..." indicates cases where transitions were not covered in the VLT/UVES spectrum. Ranges of 
values correspond to lower and upper limits measured in the cases where there was blending. 
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Table 2. Statistical Results on Very Strong and Moderately Strong Absorption Systems 





Fe IIA2374 
Mg IIA2796 


Fe 
Mg 


IIA2383 
IIA2796 


Fe 
Mg 


IIA2587 
IIA2796 


Fe 
Mg 


IIA2600 
IIA2796 


1 \\ ^ \ j t"n f~\rifM' , £i'!"ol"\7" "\rti cHvi^Tifr till 'y t l 

P(A-D) (same) 


1.7e-3 
1.6e-3 




0.4e-3 
0.05e-3 




7e-6 
0.02e-3 




2e-6 
3e-6 


P(K-S) (moderately vs very strong - z abs < 1.2) 
P(A-D) (same) 


1.3e-3 
0.7e-3 




0.2e-3 
0.15e-3 




3e-6 
0.07c-3 




0.6e-06 
4e-06 


P(K-S) (moderately vs very strong - z abs > 1.2) 
P(A-D) (same) 


0.2 
0.14 




0.014 
0.009 




0.14 
0.08 




0.07 
0.09 


P (K-S) (very strong z < 1.2 vs very strong z abs > 1.2) 
P (A-D) (same) 


0.04 
0.08 




0.015 
0.05 




0.03 
0.02 




0.005 
1.2e-3 


P (K-S) (moderately strong z aba < 1.2 vs moderately strong z aba > 1.2) 
P (A-D) (same) 


0.9 
0.6 




0.3 
0.4 




0.7 
0.6 




0.7 
0.3 


Bootstrap (moderately strong vs. very strong z aba < 1.2) 


0.04% 




0.8% 




0% 




0% 



This table provides the probability (P) that several pairs of subsamples of our data points in Figurc|6]wcre derived from the same population. 
The subsamples were created by selecting different cutoffs of 14 7 r (Mg II) (moderately strong Mg II corresponds to 0.3 < W r < 1.0 A and 
very strong Mg II absorption to W r > 1.0 A) and redshift (we divide the samples into low and high redshift using the median of the z aba 
values for each Fe II transition and taking the average of the four medians, which corresponds to a value of 1.20). We use the statistical 
tests Kolmogorov-Smirnov (K-S) and Anderson-Darling (A-D). We include the values for the four Fe II transitions analyzed (in each different 
column). Rows indicate the two subsamples in comparison, except for the last row where we include the results of a bootstrap experiment: we 
select randomly from the subsample of moderately strong absorbers at low redshift as many points as there are very strong absorbers at low 
redshift, and calculate the number of times that they cluster at high values of Wr(Fe n)/l4 7 r (Mg n) as the subsample of very strong absorbers 
at low redshift does. Whether or not they cluster was defined as whether the W r (Fe n)/Wr(Mg II) values of the randomly selected data points 
were larger than the median of the W r (Fe ll)/W r (Mg II) as many times as the values of the very strong absorbers. The value of the median 
was obtained from all the values of Wr(Fe ll)/14 7 r (Mg II) for the very strong and moderately strong subsamples combined; this value differs 
for each Fe II transition. 
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